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ABSTRACT 


We provide a detailed analysis of direct detection lightwave systems employing 
an optical preamplifier at the receiver a1 ‘ive the closed form expression for the 
bit error probability of WDM systems empioying on-off keying (OOK) as modulation 
format. In our analysis, we consider various cases in which the receiver model uses 
either a finite-time integrator or Fabry-Perot filter operating in a single channel 
or multi-channel environment. We take into account the optical amplifier noise, 


the postdetection receiver noise, the shot noise, and the effect of the nonzero laser 


linewidth. 
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I. INTRODUCTION 


Current development of optical amplifiers has advanced to the stage that their 
widespread use in lightwave systems is certain in the near future [1]. Receivers 
using optical amplifiers have been shown to be substantially more sensitive than 
their counterparts (2|. Of primary importance is the use of the optical amplifier as 
a preamplifier in a direct detection receiver. In this study, a unified approach to the 
performance evaluation of direct detection on-off keying (OOK) lightwave systems 
which take into account the effect of postdetection thermal noise, shot noise, and 
the impact of nonzero laser linewidth is provided. 

The direct detection OOK receiver to be analyzed is shown in Fig. 1. The 
received signal is amplified by an optical amplifier of gain G. The optical amplifier 
introduces additive white Gaussian noise (AWGN) with zero mean and power spec- 
tral density (PSD) No/2 where No = N,,hf(G — 1) [2, 12, 33]. The parameter Nz, 
represents the spontaneous emission factor which is unity for an ideal amplifier, A 
is Planck’s constant (6.626 x 10-*4 J.s), and f is the frequency. Because the optical 
amplifier is not polarization independent, a polarizer is needed to pass the desired 
signal and to block the light in the orthogonal polarization. In this study, the optical 
bandpass filter is modeled by both a finite-time bandpass integrator with integration 
time T’ [6, 9-11, 33] and a Fabry-Perot filter. 

For the finite-time integrator, the noise bandwidth 1/T”’ of the bandpass in- 
tegrator is chosen to be the same as the noise bandwidth of the optical filter. The 
equivalent lowpass impulse response of the integrator is h(t) = (1/T’)[u(t)—u(t-T’)| 
where u(t) is the unit step function. For convenience of analysis, it is assumed that 


the bit duration T is a multiple integer of T’, that is, T= MT’ [33]. 
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OOK receiver structure. 


Figure 1: 


The photodiode has a responsivity R = neq/hf where 7. < 1 is the quantum 
efficiency and q is the electron charge (1.6 x 10-!9 C). The output current of the 
photodiode, which is proportional to the squared envelope of the receiver signal, is 
processed by a detector which consists of a low noise amplifier, an integrate-and- 
dump filter and a slicer (33). 

In the case where a finite time integrator is employed, the amplifier noise 
effect is most pronounced at low postdetection thermal noise. With a higher gain, 
the amplifier noise is dominant and, consequently, the net amplifier gain is reduced 
considerably. In the case where a Fabry-Perot filter is employed as an optical filter 
at the receiver model, the equivalent lowpass impulse response of the Fabry-Perot 
filter can be well approximated by an RC filter for both single ar? multi-channel 
within the frequency range |f — fo| < F.SR/207 where FSR is the free spectral 
range and fo is the center frequency of the Fabry-Perot filter. For example, given 
FSR = 3800 GHz, the approximation works very well for |f — fo| < 60.5 GHz. In 
other words, the effects of adjacent channels within 121 GHz can be included. The 
same model without the optical amplifier is analyzed in [29] and it is seen that this 
model enables use to obtain a closed form analytical expression for the bit error 
probability. In all bit error probability derivations, whether a finite-time integrator 
or a Fabry-Perot filter is employed as the optical filter in the receiver model, it is 
assumed that all channels are bit synchronous as in (29, 30]. 

This thesis is organized into three sections. In the second section, the receiver 
models employing finite-time integrator, and the analysis is a reproduction of [33], 
Fabry-Perot filter for single channel and Fabry-Perot filter for multi-channel are 
analyzed. Each of these analyses obtain the mathematical framework, bit error 
probability, and the numerical results deriving the detected signal erivelope and its 


density function. The last section provides the conclusions. 
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II. ANALYSIS 


The analysis is organized into three subsections. All three provide the math- 
ematical framework, bit error probability derived from the statistics of the white 
Gaussian postdetection thermal noise and the numerical results for the receiver 


model considered. The receiver model to be analyzed is shown in Fig. 1. 


A. FINITE-TIME INTEGRATOR MODEL 
The optical bandpass filter is modeled as a finite time bandpass integrator 
with integration time T’. The noise bandwidth 1/T’ of the bandpass integrator is 
chosen to be the same as the noise bandwidth of the optical filter. This model is 
used because it is analytically tractable. This section is taken from (33]. 
1. Mathematical Framework 
For mathematical convenience we adopt the complex envelope notation 
of a real OOK signal. Thus, for a given transmitted bit 6; the corresponding received 


signal at the input of the optical bandpass filter in Fig. 1 is designated as follows: 
r;(t) = VGP bye? pp(t) + n-(t) + jne(t) (1) 


where pr = 1,0 <t < T and zero otherwise, G is the optical amplifier gain, P is the 
peak power at the input of the optical amplifier, 6(t) is the OOK laser phase noise, 
and n,(t) and n,(t) are the independent in-phase and quadrature components of 
the additive zero mean white Gaussian noise representing the amplifier spontaneous 
emission noise. The PSD of n,(t) and n,(t) is No/2 (we use No/2 instead of the usual 
lowpass PSD No because the magnitude of the signal component is /GP b; instead 
of the usual /2GP }; for complex envelope). Let r be the extinction ratio of the 


laser source defined as the ratio of the transmitted power for the logical zero to that 


5 


for the logical one. Then bg = \/r/(1 +1) for the logical zero, and 6, = yi/(l+r) 


for the logical one. 


a. Approximation of the Detected Signal Envelope 


The output signal of the optical bandpass filter is given by 
1 se . 
s(t) = Ti fut) — u(t —T'\ri(t — r)dr 


1 ft 
= xf rj(t — r)dr 
1 ft d 
= af rinar (2) 
The photodetector with responsivity R detects the squared envelope Rjs;(t)|?. This 
signal current plus the shot noise current v;(¢) generated by the photodiode, and 


the postdetection thermal noise current w(t) are scaled by 1/TRGP and integrated 


to produce the decision variable 


Vee Pera [ rudoracs ate [vunat 
= X;+V+w (3) 
here 
ee sap [Renee (4) 
Yo= ROP [ena (5) 
W = TRGP [ w(eyat (6) 


We model w(t) as a zero mean white Gaussian noise current of spectral density Wo. 
It is seen that the variance of o?, of w(t) is 


Wo 
T R?G? Pp? (7) 


oy, = 


6 


The shot noise current v,(t) can be modeled as a zero mean Gaussian noise current 


with spectral density Vo; as follows [15]: 
Voi = gRE{|s;(t)|?} (8) 


The variance of, of v,;(t) is thus given by 


oF, = TRoeprE A lsi(t)?} (9) 


The value of of, can be computed from (A1) and (A4) as follows: 


0, = —aMi _ (FBT, norm _ 1) , IMNewAf(G —1) 
“ ~ TRGP(r AT)? T? RG? P? 


(10) 
where £ is the laser linewidth. 

The statistics of X; are very difficult to obtain. Therefore, we make 
the following approximation of X,;. We note that the equivalent lowpass impulse 
response h(t) of the optical bandpass filter is assumed to be time-limited to T’ = 
T/M. Thus, the integrand R\s;(kT’)|?, k = 1,2,---, M in (4) is evaluated at t = kT’ 
where there are & independent and identically distributed (iid) samples of R|s,(t)|? 
during an interval T. Thus |s;(k7’)|? is obtained from disjoint integration intervals 
where the phase noise @(t) has independent increments and the noise n,(t) and n,(t) 


are altogether independent over disjoint intervals. Using the above fact we can 


approximate X; as the sum of M iid samples T’R|s;(kT’)|?/TRGP. That is 


Wy 42 
Xi Tere 3 |si(kT’)| (11) 
Substituting (1)-(2) into (11) we obtain 
2 
X; & 5 |e Ae Nan + 5 (12) 
=1 
where 
~ Ls ia) 
Lh = [ ey Oat (13) 
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l kT" 


Nx Oe -(t)dt 14 
: T'/MGP (evr ( ) ( ) 
N, ee eee ss d 15 
a 3 t)at 
, T'/MGP | ( ) ( ) 


The random variables N,, and N,,; are independent Gaussian random variables with 


zero mean and variance o? given by 


ya No _ Nohf(G-1) 
2TGP 2TGP 


(16) 


The approximation of X; in (12) is equivalent to the modeling of the integrate-and- 
dump filter as a discrete-time integration that sums over M samples taken every T’ 
seconds at the output of the photodetector [6, 13]. We remark that both X; and its 
approximation have the same mean value. If |s;(t)|? has a constant spectral density, 
then X; and its approximation also have the same variance. This happens when M 
is large and the spectral density of |s;(t)|? can be considered constant within the 
bandwidth 1/T of the integrate-and-dump filter. 
b. Probability Density Function of the Detected Signal 

Envelope 

From (11) we observe that X; is approximated by the sum of squares 
of 2M iid Gaussian random variables. Therefore, X; is approximately noncentral 


chi-square distributed with the following conditional probability density function 
(pdf) fx;(zily) [18] 


1 ee VANE x, + b?y/M 
fx;(zily) = 202 (zn) exp {=e 


a 


TGP (Mere { -POPtat + ait) 
Naphf(G—1) \ by PL NB AFG - 1) 


Iu-1 
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F 2TGP,/byz;/M 7 
M-1\ “NhF(G —1) a) 
where Iy_,(-) is the modified Bessel function of order M — 1, and 7 is the value 
assumed by the random variable defined as follows: 
M 
P= ylZ/ (18) 
k=1 
In the special case when the extinction ratio r = 0, then bp = 0, and fx, (zo) in (17) 


is a chi-square probability density function (pdf) given by [18] 


1 TGP. Vp TGP 2, 
Fxo(20) = Ty (acer) 0: oo | - RED (19) 


The pdf of fx,(zi|7) can be obtained by averaging fx,(z;|7) over 7. Thus knowing 
fr(7) we get 
fx (ai) = [ faslzily) Solver (20) 
2. Bit Error Probability 
For a detection threshold a, the bit error probability conditional on the 
mean Zo and xz, assuming a combined shot noise current and postdetection thermal 
noise current spectral density Vo; + Wo is given by 


P,(z09, 21) = zerfe (52) + zerfe (-*) (21) 


1 


where the complimentary error function erfc(-) is defined as 


7 hee 
erfe(a) = al e~* dr (22) 


and o}., ? = 0,1 are the variances of Y; in (3) 


oY, = o3, + ow 
pees ee ayy 4 —We _ 


The bit error probability P, is obtained by taking the expectation of P,(zo,2,) with 
respect to ro and 2, by using (17) 


P= a erfc (5 fx,(xo)dzo 


+h erfe (25-*) eit (24) 


The optimal threshold that minimizes the bit error probability is the value that 


satisfies the following equation 


00 a-2 20 m—-2 
[ erfe(5*2) fxo(ao)de0 = [ erfe (==) Fx, (21)dz, (25) 


From the above analysis it is seen that the evaluation of the bit error probability P, 


requires the knowledge of the pdf fr(7) of I in (18) which in turn requires the pdf 
of |Z;,|? in (13). The pdf of the random variable |Z,|? has been studied extensively 
in (13, 14], [20, 21]. The evaluation of P, in (24) is computer intensive even in the 
ideal case of no laser phase noise. When laser phase noise is taken into account, the 
pdf of must be obtained from the pdf of |Z,|* via an M-fold convolution. 
3. Numerical Results 

In this section we present numerical results for a direct detection light- 
wave system with the following parameters: r = 0.05, 1/T = 500 Mb/s, G = 100 
and 1000, and Wo = 10-™%, 10-*3, 10-2? (A?/Hz). Figure 2 shows the bit error 
probability P, versus the input power P (dBW) as a function of postdetection noise 
spectral density Wo for a receiver without the filter and the optical amplifier. Figure 
3 shows P, versus P for G = 100, M = 10 without phase noise and with phase noise 
of ST = 1 as a function of Wo. At P, = 10-", the net amplifier gain without (with) 
phase noise is 19.6 dB (19 dB) for Wo = 10-?? A?/Hz, 18.6 (18 dB) for Wo = 10-* 
A?/Hz, and 16.3 dB (15.7 dB) for Wo = 10-4 A?/Hz. It is seen that the amplifier 


noise effect is most pronounced at low postdetection noise. 
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Figure 2: Bit error probability versus input power for a system without optical 
amplifier as a function of postdetection noise spectral density Wo. A: Wo = 107-72 
A?/Hz, B: Wo = 10-* A?/Hz, C: Wo = 107% A?/Hz. 
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Figure 3: Bit error probability versus input power as a function of postdetection 
noise spectral density Wy with G = 100, M = 10. A: Wo = 10-* A?/Hz, B: 
Wo = 10-* A?/Hz, C: Wo = 10-* A?/Hz. 


12 


Figure 4 shows P, versus P for G = 1000, M = 10 without phase noise 
as a function of Wo. At P, = 10-", the net amplifier gain is 26.4 dB at Wo = 107” 
A?/Hz, 22.1 dB at Wo = 10-”? A?/Hz, and 17.2 dB at Wo = 107*4 A?/Hz. It is 
seen that with a higher gain, the amplifier noise is dominant and consequently, the 
net amplifier gain is reduced considerably. Comparing the two systems that employ 
optical amplifiers with gain G = 100, and G = 1000, respectively, we observe that 
i is a net improvement of 6.8 dB at Wo = 10-2? A?/Hz, 3.5 dB at Wo = 10- 
ro and only 0.9 dB at Wo = 10-*4 A?/Hz for G = 1000 over G = 100. Besides 
that, Fig. 5 shows the optimized threshold versus input power P for G = 100, 
M = 10, 30, 50, and 100. As it is described in (25), the optimal threshold that 
minimizes the bit error probability is the value that satisfies the equation given 
in (25). As we can see from Fig. 5, optimized threshold value decreases as input 
power increases for all M values. For considered values of M, which are 10, 30, 50, 
and 100, the optimized threshold appears to converge to about 0.5 for higher input 
power values. We conclude that a larger gain amplifier should be used when the 
postdetection noise is large, and a smaller gain amplifier should be used when the 


postdetection noise is small. 


B. FABRY-PEROT FILTER FOR SINGLE CHANNEL 

The desired signal is amplified by a preamplifier and then filtered by using a 
Fabry-Perot filter. The photodetector has a responsivity R (A/W). The detected 
current is amplified by a low noise amplifier that adds a postdetection thermal 
noise with spectral density No (A?/Hz). The decision variable at the output of the 


integration is compared to the threshold in order to be able to determine whether 


a bit zero or bit one is present. 


Bit Error Probability 
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Figure 4: Bit error probability versus input power as a function of postdetection 
noise spectral density Wo with G = 1000, M = 10. A: Wo = 10-*? A?/Hz, B: 
Wo = 10-* A?/Hz, C: Wo = 107% A?/Hz. 
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Figure 5: Optimized threshold versus input power as a function of M with G = 
1000, postdetection thermal noise Wp = 10-2? A?/Hz. 
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1. Mathematical Framework 

Again, we adopt the complex envelope (equivalent lowpass) notation of 
a real OOK signal for mathematical convenience. For single channel analysis, we 
normally ignore the intersymbol interference (ISI) and adjacent channel interference 
(ACI) components. Detailed analysis including ISI and ACI components are per- 
formed in the multi-channel case. For a given transmitted bit b;, the corresponding 
received signal at the input of the Fabry-Perot filter in Fig. 1 is designated as 
follows: 


ri(t) = VGP boopr(t) + ne(t) + jna(t) (26) 


where G is the optical amplifier gain, P is the peak power at the input of the 
optical amplifier, n.(t) and n,(t) are the independent in-phase and quadrature com- 


ponents of the zero mean AWGN representing the amplifier spontaneous emission 


noise (ASE) and 
1 O<t<T 
Pr = (27) 


0 otherwise 


The Fabry-Perot filter can be characterized by the following equivalent lowpass 


transfer function [31, 32] 
_ l-p 1-—A-p 
Hf). = 1 — pe-27/FSR 1 — p 
(28) 


H(f) = 7s = 


1 — pcos (75) + jpsin (75) cages 
where p is the power reflectivity, A is the power absorption loss (A = 0 for an ideal 


Fabry-Perot filter) and FSR is the free spectral range. For |f| < FSR/20 and 
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assuming A = 0, we can approximate H(f) as follows 


l-p 1 
Af) A= 
2a fp ._anfp 
Q-)+izce 1+Iq—))FSR 
(29a) 
1 
——FxF ? lf| < FSR/207 
Lj 
c 
where 
oe (29b) 


The free spectral range FSR can be related to the full width at half 

maximum (FWHM) bandwidth B and the finesse F of the Fabry-Perot filter as 
rsp= "VP". BP (30) 
Thus, if the signal is bandlimited to |f| < FSR/207, we can actually use (29), but 
with the frequency covering the entire frequency spectral range, that is, we can truly 
approximate (28) with a single-pole RC filter with the following tra::sfer function 


and impulse response, respectively, 
H(f) = ——ey (31) 


h(t) = ce*, t>0 (32) 


The magnitudes of H(f) of the Fabry-Perot filter in (28) and its approximated single- 
pole RC filter in (31) for p = 0.99, F = 312.6, B = 12.16 GHz, and FSR = 3800 
GHz, remain identical and attenuate rapidly while the phases differ as the frequency 
increases [29]. 

The above approximation is valid for dense wavelength division multi- 
plexing (WDM) analysis when the filter finesse F is large or equivalently the FWHM 


bandwidth B is small since the equivalent lowpass signal must be bandlimited to 
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about |f| < FSR/207. When |f| > FSAR/207, the magnitude of H(f) is very small 
and therefore, the effect of ACI beyond this frequency range is negligible. 


a. Derivation of the Detected Signal Envelope 
Since we are interested in the detected bit bo in the time interval 
(0,7), we consider the output preamplified and filtered S;(t), given by 
S,(t) = Sp(t) 0<t<T (33) 


where S(t) is the desired signal, which is the convolution of the desired signal and 


impulse response A(t) in (32), as 


h(t) = ce t>0 (34) 
Sa(t) = VGP boo [we-nar (35) 
Sp(t) = VGPboo(l-e*), O<t<T (36) 


The photodetector with responsivity R detects the squared envelope R]S;,(t)|?. This 
signal plus the postdetection thermal noise current w(t) are scaled by 1/TRGP and 


integrated to produce decision variable, 


~ 2 fpf pista 4 —— |" 
% = ppop [ RISi(t)Pdt + sas [ w(t)dt (37) 
Y= X;+W (38) 
where 
1 T 2 
Xi = apap f Rlsorat (39) 
1 T 
W = zone | w(t)dt (40) 


We model w(t) as a zero mean AWGN current of spectral density of Wo. It is seen 


that the variance of 07, of w(t) is 


WoT Wo 


*W = TaURGPP ~ TIRGPF o 
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As it is described in (37), the decision variable Y; can be easily found for zero and 


one bit as follows, 
Yo Xo+W (42) 
NY Mi+W (43) 


As it can be seen from above, since we ignore shot noise in this analysis, the variances 


of Y; and Y2 are equal to each other. 
OF, = oY, = (44) 
X; can be computed from (B.5) and found as (B.6) 
Xi = Wo [1 — eer ee aa) (45) 
ae cT 2cT 


Using the fact that the noise n,(t) and n,(t) are altogether independent over disjoint 
intervals and the random variables N, and N, are independent Gaussian random 
variables with zero mean and variance a? computed from (B.18), we can define 


cy 


variance as 
N,yhs(G —1)c 
2 _ i‘spi'f 
“1 = "4MGP 46) 
b. Derivation of Probability Density Function of the Detected 
Signal Envelope 
As it is described in the first section, X; is approximately noncentral 


chi-squared distributed with the following probability density function (pdf) f,;(zi) 
[18], 


M-1/2 


oF 


y2 


2 1 
: er 2 _ 211 — 7 _— p-2cr 
die Po [: pl ey aT 7 ) 
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ae for bit one (49) 


r a 
Vien for bit zero (50) 


h;(G = 1)e 
4MGP (51) 


In the special case where the extinction ratio r = 0, since boo = 0 for bit zero, pdf 


of f,,(z;) in (47) becomes a chi-square pdf given by [18], 


fio(20) = onan” _~ 0/203 (52) 


2. Bit Error Probability 
For a detection threshold a, the bit error probability is conditional on 


the mena zo and z, and given by 


7 1 a-—2% 1 %1-a 
P,(t0, 71) = gorse (==) + 7 fc (a=* =) (53) 


where erfc(-) is defined as in (22) and o},, i = 0,1, are the variances of Y; in (37) 
OY, = OY, = OF, = Ow (54) 


We can represent both variances as o?, since they are equal to each other. The bit 
error probability is obtained by taking the expected value of P,(zo, 21) with respect 


to Zp and 2, using (47), 


R=; -f erfe( Se +2) fay(to)dt0 + = f° erfe (Hes ae =) fe(ti)dz, (55) 


The optimal threshold a that minimizes the bit error probability is the value that 


satisfies the following equation also given in (25) 
I erfc Ga *) fzodzo = i erfc (a =) fe, (21)dz, (56) 
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3. Numerical Results 

In this section, we present the numerical results for a single channel 
direct detection lightwave system employing a Fabry-Perot filter and an optical 
amplifier at the receiver model in Fig. 1 with the following parameters: bit rate 
R, = 500 Mbps, optical amplifier gain G = 100 and 1000, photodiode responsivity 
R=0.5, FSR = 3800 GHz, c = 38.4 GHz, p = 0.99, thermal noise spectral density 
Wo = 107??, 10725, 10-*4 (A?/Hz), and M = 10, 13, 20, 40 where M = T/T’, 
T is the bit duration (1/R,) and T’ is the time constant of the optical bandpass 
filter. Since bit rate R, = 500 Mbps, bit duration T is 2 ns. If we chop the impulse 
response of the filter at 0.05 ns, 0.1 ns, 0.15 ns, and 0.2 ns, M has the values of 
40, 20, 13, 10, respectively. Figures 6, 7, and 8 show the bit error probability P, 
versus peak input power P (dBW) as a function of M for the thermal noise spectral 
density Wo = 107??, 10-25, 10-74 A?/Hz with G = 100. As M decreases, it is 
obvious that more input power is required to achieve P, = 107 as compared to 
finite-time integrator. In the case where finite-time integrator is employed as M 


decreases, less input power is required because the noise variance given in (16), 


3? = Neh f(G-1) 
2GPT 


is independent of M and even if T’ is varied, it has constant T. Therefore, varying M 
does not effect the noise variance but it effects the signal power. For the Fabry-Perot 


filter, the variance given in (46), 


2 Nwphf(G— Ve 
°: =~ 4MGP 


depends on the M value, and as T’ is varied, M varies accordingly (M = T/T’). 
So, as M decreases, noise variance gets larger and more power is required to achieve 


P, = 107. As it can be easily seen from these figures, the role of T’ is more 
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Figure 6: Bit error probability versus input power as a function of M with G = 
100, Wo = 10-”* A?/Hz (Fabry-Perot filter). 
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Figure 7: Bit error probability versus input power as a function of M with G = 
100, Wo = 10-8 A?/Hz (Fabry-Perot filter). 
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Figure 8: Bit error probability versus in 
100, Wo = 10-4 A?/Hz (Fabry-Perot filter 
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a power as a function of M with G = 


significant for smaller thermal noise values. This is more pronounced with a high 


gain amplifier (G = 1000) as can be seen in Tables 1 and 2 and Figs. 9, 10, and 11. 


Detailed comparison is as follows: 


able 3) Comparison between M = 10 and M = 40 with G = 100 (see Figs. 6, 
, and 8). 


Table 2: Comparison between M = 10 and M = 40 with G = 1000 (see Figs. 9, 
10, and 11). 


—78.7 dB —73.9 dB 4.8 dB 


G = 1000 


10 


Bit Error Probability 


10° M = 10 


-82 -81 -80 -79 -78 -77 -76 -75 -74 -73 -72 
Input Power (dBW) 


Figure 9: Bit error probability versus input power as a function of M with G = 
1000, Wo = 10-2? A?/Hz (Fabry-Perot filter). 
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Figure 10: Bit error probability versus input power as a function of M with 
G = 1000, Wo = 10-8 A*/Hz (Fabry-Perot filter). 
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Figure 11: Bit error probability versus input power as a function of M with 
G = 1000, Wo = 10-4 A?/Hz (Fabry-Perot filter). 


When we compare the finite-time integrator and Fabry-Perot filter with 
each other, Table 3 shows that the finite-time integrator has better performance than 
the Fabry-Perot filter for larger T’. This improvement of the finite-time integrator 
over the Fabry-Perot filter gets more significant with high gain amplifiers and lower 
thermal noise values as displayed in Table 4. Performance margins decrease with 
increasing thermal noise values. 

Figures 14 and 15 show that the Fabry-Perot filter has better perfor- 
mance than finite-time integrators with a low amplifier gain (G = 100) for low 
thermal noise values (Wo = 10-*% A?/Hz, 10~** A?/Hz) and both have the same 
performance for a higher thermal noise value of Wo = 10-?? A?/Hz with smaller T’, 
which implies a higher values of M due to the relation of M and T’ (M = T/T’) for 
a constant bit rate. This improvement of the Fabry-Perot filter over the finite-time 
integrator is more pronounced with a high gain amplifier (G = 1000) as shown in 
Figs. 16 and 17 and tabulated in Tables 5 and 6. 

As a result, the choice of T’ and consequently, the value of M plays 
a great role for the systems employing the Fabry-Perot filter and the high gain 
amplifier when the postdetection thermal noise is small. The systems employing 
finite-time integrators and high gain amplifiers have better performance than those 
employing the Fabry-Perot filter when 7” (time constant of the filter) is large and 
thermal noise is small. The systems with the Fabry-Perot filter and the high gain 
amplifier have better performance than the finite-time integrator when T’ is small 


and thermal noise is low. 
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Table 3: Comparison between FTI and FP for G = 100, r = 0, and M = 10 (see 
Figs. 6, 7, 8, and 12). 


) 
Remarks 
Improvement of 
Wo (A? / Hz) 


Table 4: Comparison between FTI and FP for G = 1000, r = 0, and M = 10 (see 
Figs. 9, 10, 11, and 13). 
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Figure 12: Bit error probability versus input power as a function of postdetection 
thermal noise Wo with optical amplifier gain G = 100, extinction ration r = 0, 
M = 10. A: 10-”* A?/Hz, B: 10773 A?/Hz, C: 107*4 A?/Hz (finite-time integrator). 
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Figure 13: Bit error probability versus input power as a function of postdetection 
thermal noise Wo with optical amplifier gain G = 1000, extinction ration r = 0, 
M = 10. A: 10-” A?/Hz, B: 10-% A?/Hz, C: 10-4 A?/Hz (finite-time integrator). 
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Figure 14: Bit error probability versus input power as a function of postdetection 
thermal noise Wo with optical amplifier gain G = 100, extinction ration r = 0, 
M = 50. A: 10-?* A?/Hz, B: 10-7 A?/Hz, C: 10-74 A?/Hz (finite-time integrator). 
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Figure 15: Bit error probability versus input power as a function of postdetection 


thermal noise Wo 
M =50. A: 107? 


with optical amplifier gain G = 100, extinction ration r = 0, 
A?/Hz, B: 10-*° A?/Hz, C: 10-4 A?/Hz (Fabry-Perot filter). 
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Figure 16: Bit error probability versus input power as a function of postdetection 
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thermal noise Wo with optical amplifier gain G = 1000, extinction ration r = 0, 
M = 50. A: 107%* A?/Hz, B: 10-73 A?/Hz, C: 10774 A?/Hz (finite-time integrator). 
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Figure 17: Bit error probability versus input power as a function of postdetection 
thermal noise Wo with optical amplifier gain G = 1000, extinction ration r = 0, 
M = 50. A: 10-” A?/Hz, B: 10-79 A?/Hz, C: 10-74 A?/Hz (Fabry-Perot filter). 
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Table 5: Comparison between FTI and FP for G = 100, r = 0, and M = 50 (see 
Figs. 14 and 15). 


Table 6: Comparison between FTI and FP for G = 1000, r = 0, and M = 50 (see 
Figs. 16 and 17). 


Remarks 
Improvement of 
Fabry-Perot is 
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C. FABRY-PEROT FILTER FOR MULTI-CHANNEL 


In the previous section, we considered only the desired bit signal rather than 
the ISI and ACI components since we dealt with a single channel. In this section, 
we take into account the ISI and ACI components and our analysis focuses on the 
multi-channel. 

1. Mathematical Framework 
For convenience, we designate channel 0 as the desired channel and chan- 
nel k as the adjacent channel where k = —M/2,---,—1,1,---,M/2 and M is an 
even integer. We consider the equivalent lowpass (complex envelope) data signal in 


channel 0 and channel k as follows, 


0 

to(t) = 30 boapr(t iT) (57) 
t=—bo 

b(t) = 3 beee™pr(t — £7) (58) 
i=-L 


where T is the bit duration, bo;¢{0,1} is the bit in channel 0 in the time interval 
(7, (t+1)T), by ce{0, e7*} is the 2 bit in channel k in the time interval (¢T, (€+1)T), 
¢, is a phase offset between channel 0 and channel k, which is assumed to be 
uniformly distributed in (0,27) radians, and w, is the frequency spacing between 
channel 0 and channel k in radians. 
The function pr(t — 27) is defined as 
1 iT<t<(4+1)T 
pr(t~ iT) = (59) 


0 ~— otherwise 


The non-negative integers Lo and L in (57) and (58) represent the number of bits 


in channel 0 and k, respectively, that precede the detected bits boo. The received 
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WDM equivalent lowpass signal at the input of the Fabry-Perot filter is given by, 


M/2 
r(t)=VGPb&(t)+ >. VGP&(t) (60) 
ka M/2 


kyo 
where P is the received optical power and G is the gain of the preamplifier. 

As we discussed in the previous section, the impulse response of the 
Fabry-Perot filter is approximated by a single-pole RC filter of which the impulse 


response is given in (32) as, 
h(t) = ce~, t>0 (61) 


a. Derivation of the Detected Signal Envelope 


The output preamplified and filtered signal S;(t) is given as, 
Si(t) = Sp(t) + Srs1(t) + Saci(t), 0<t<T (62) 


where Sg(t) is the desired bit signal, $;5;(t) is the intersymbol interference signal, 
and S4c;(t) is the adjacent channel interference signal. These signals are evaluated 


to obtain output filtered signals by using (60) and (61), 


Set) = VGP boo [ "alt — 1) de 
= VGPho(l—-e“), O<t<T (63) 
—t i+ 
Stsi(t) = VGP Jo bo, de FRG ae 


t=Lo 


-1 
VGPe~* > by; (eT eT), Ot ST (64) 
t=Lo 


-1 
YS bee | ie h(t — nena] 


= — ZL rT 


M/2 
Sucilt) = VGP 5 | 


ku—M/2 


kd 


t ‘ 
the | h(t — near 


c+ju, =-L 


M/ ~ 
= VGPe* 3 ! | > bie (eletimaherur = cosine 
é 


hu —M/2 
kg 


+bp,0 (e(+7*) — 1) } , O<t<T (65) 


The Fabry-Perot filtered output S;(t) is detected by the photode- 
tector. The detected current signal is amplified by a low noise amplifier (LAN) that 
contributes a postdetection thermal noise w(t) with spectral density No (A?/Hz). 
The decision variable at the output of the integration is compared to a threshold a 
to determine whether bit zero or bit one is present. The decision variable Y is given 


by (37) as, 


Y=X;+w (66) 
. where 
— 1 a . 2 
Xi = apap if RisePat (67) 
1 T 
W = apap ‘| W(t)dt (68) 


We obtain the signal component X; by substituting (62), (63), (64), and (65) into 
(67) as a function of cT and w,.7, which represents the impact of ISI and ACI, 
respectively. X; can be obtained from (C.19) for the worst case analysis as follows 


[29], 


2 1 
. = 2 —cT7 —2cT 


1 
TRGP 


40 


+52 (a (jo e*T) + bogb- (1 — 26-7 + e°¥7)| 


+ [ioe + Za — er) Re(0}(6. — bo) 


= Ty) (69a) 


ky 
where 


boo = 0 for bit zero boo = 1 for bit one, 
b. =1 for bit zero b. =0 for bit one, 
b=1 for bit zero b=0 for bit one 
in the worst case analysis. 
For the exact case analysis, X;, which is derived in [29] for the case 


without the optical amplifier, can be modified for our case as follows, 


X; = saa | RIS;,(t)[Pdt 
2 _d 1 By 
Xi = aaa [PRGPHo[1- a (1-e 7) +55 (1-e**)| 


2 
+TRGPS = (1 = gt) Ee j (eter — a) 


+TRGP (1 - e-*7) 


bee 
3 & 0+ jw,T /cT 
M/2 


2 
(erent Hee) 


1 M2 MR by.05%, 0 
+TRGP- > (1+ juxT/cT)(1 — jwmT/eT) 


ku~M/20 ma-M/2 
kyo mx 


4] 


—(cT-jw,T) _ —(cT+jwmT) _ 
cT Sf ie) 
1 — ju,T /cT l+jwaT/cT 2 


errs LR. M/2 M/2 —1 bh 0° 0 
2 eum 
7 > x 25 (1 + jupT/cT)(1 — jwmT/cT) 


i -M/2) ms—-M/2 


kyo m0 


(elte entry? = eee) 1 — e(Ttiue7) = 1 (1 a enn) 
l+juwaT/cT 2 


+GPT boo (1 +e 2eT _ 2e~*T) > bo (etrner = eT) 
i=-Lo 


1 —2cT -c 
$GPT ino ZR (1 + e7*T _ 2e~“7) 


Mie: ak b ; 
5 > ke (elcThiveT (C41) mt eleTrieaT)) 
a 1 + jw.T /cT 


ku —M/2 
kgd 
42 > bio et eel gent 1 _ 1 ter 
1+ jw.T /cT er 2 
ku—-M/2 
kx 
e-(cT-juaT) _ 1 GPT 
ee (i+1)cT _ ,icT 
= 1 — ju, T /cT )h+ oT Pay (c af ) 


R. on > ad See (1 - oe) 
1+ jw,T/cT 


t=-L 


feletssmaT Ker?) ot e(cT+iu. T)¢) 


M/2 b (1 - —(cT-ju,T) 
k,0 e ) —2T 
+ 25. Ts ne ( i=jurie  S i) x52) 
ku —M/2 
kd 


Using the fact that the noise n,(¢) and n,(t) are independent over disjoint intervals 


and the random variables N, and N, are independent identically distributed (iid) 
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Gaussian random variables with zero mean and the variance o/ computed from 
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(B.18) and given by (46), we can employ the variance as, 


32 = Neohs(G = 1c 


; TGE as described in (46) (70) 


b. Derivation of Probability Density Function 
As it is explained in the first and second section as well, X; is 


noncentral chi-square distributed with the following pdf also given in (47) 


(\ M-1/2 , 
fx,(Xi) = sa (4) e7 (it ¥*)/207 hes (=) (71) 


where 1? is described as in (69) 


1 : 
boo = (+r) for bit one (72) 
bo = aa for bit zero (73) 
Naphs(G — 1)c ‘ 
an pif 
o% = —Guap (46) (74) 


The special case where the extinction ratio r = 0 introduces the chi-square pdf given 
by (52) and [18], 


fieo(20) = awe ~~ 0/208 (75) 


2. Bit Error Probability 
For a detection threshold a and an ISI/ACI pattern b = {bo;, 5, .} where 
t= —Lo,---,0 and k = —M/2,---,M/2(k # 0), the bit error probability is defined 
as in (53), 


_1 a — Zo(b) 1 r,(b)—a 
Aamo (Ses) tae (ae) sf 
where 
2 Wo 
oY = ow = TTRGP) on 
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Bit error probability for the worst case analysis and exact analysis are obtained by 
taking the expected value of P,(zo,z;) with respect to zo and 2, by using (69a) and 
(69b), respectively. 


P, for the worst case analysis: 


Bad [O erte (A) steordeo + 3 [ete (BA=*) saerydnr (18) 


P, for the exact analysis: 


Fisieseie ui [et Kod zolb, ak mane 


The optimal threshold a that minimizes the bit error probability is the value that 


satisfies the following equation, 


3. Numerical Results 

In this section, we present the numerical results for the system employing 
the optical amplifier and the Fabry-Perot filter in the receiver model for multi- 
channel with the following parameters: bit rate R, = 2.56 Gbps, optical amplifier 
gain G = 1000, free spectral range FSR = 3800 GHz, M = 3 and 4, channel 
spacing J = 8, 12, and c = 38.4 GHz. In our analysis, we consider that the signal is 
band-limited and incorporate the degradation caused by the signals in the nearest 
adjacent four channels. Since we deal with multi-channel in this section, we will see 
the effect of ISI and ACI components on the performance of the system, unlike the 


single channel case. 
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For the worst case analysis with optimum threshold, boo = 0 for bit 


zero, by9 = 1 for bit one. Figure 18 shows the effect of ISI and ACI for M = 3, 
cT = 15, and J = 8, ‘2 in terms of bit error probability and input power (dBW). As 
can be seen from the graph, the single channel requires —68.8 dBW in order to be 
able to achieve bit error probability P, = 10-'®. The system with ISI but no ACI, 
requires —67.5 dBW to achieve the same bit error probability, which corresponds 
to a 1.3 dBW power penalty as compared to the single channel case. When the 
adjacent channels are taken into account for channel spacing J] = 12 and 8, the 
power penalties are 3.8 and 6.0 dBW, respectively. As the channel spacing between 
channels gets larger, the effect of ACI becomes less significant. 

Figures 19 and 20 show the effect of the value of M, which is the ratio 
of bit duration (T) to the time constant of the filter impulse response (J’). As 


explained in section 2, the noise variance has a dependency on M as given in (46), 


52 — Neh f(G=1) 
’ 4MGP 


A M decreases, the noise variance gets larger and more noise power is obtained 
accordingly. As can be seen from the figures, more input power is required for 
M = 3 rather than M = 4 at P, = 107) as expected from the results obtained in 
the single channel case in section 2. 

When we compare the worst case analysis and exact analysis with each 
other, as we expected, exact analysis requires slightly less power than worst case 
analysis in order to obtain the same probability of bit error P, = 10-!° as shown in 


Fig. 21. 
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Figure 18: Bit error probability versus input power for Fabry-Perot filter with 
G = 1000, M = 3, Wo = 10-!? A?/Hz (multichannel worst case analysis). 
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Figure 19: Bit error probability versus input power for a Fabry-Perot filter as a 
function of M with G = 1000, Wo = 10-** A?/Hz, I = 8. 
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Figure 20: Bit error probability versus input power for a Fabry-Perot filter as a 
function of M with G = 1000, Wo = 10-7? A?/Hz and J = 12. 
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Figure 21: Bit error probability versus input power for Fabry-Perot filter with 
G = 1000, Wo = 10-5 A?/Hz, M = 3, J = 8 (comparison between exact and worst 
case analysis for multichannel). 
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III. CONCLUSIONS 


We have presented a detailed analysis of direct detection lightwave systems 
employing an optical preamplifier at the receiver and derived the closed form ex- 
pression for the bit error probability of WDM systems employing on-off (OOK) as 
a modulation format. We have considered various cases in which the receiver model 
uses either a finite-time integrator or Fabry-Perot filter operating in a single chan- 
nel or multi-channel environment. We have taken into account the optical amplifier 
noise, the postdetection receiver noise, the shot noise, and the effect of the nonzero 
laser linewidth. 

For the finite-time integrator, we conclude that the amplifier noise effect is 
most pronounced at low postdetection thermal noise. A larger gain amplifier should 
be used when the postdetection thermal noise is dominant, while a smaller gain 
amplifier should be used when the postdetection thermal noise is small and amplifier 
noise dominates. 

The performance of the system employing the Fabry-Perot filter in a single 
channel environment depends on the value of M, which is defined as the ratio of 
bit duration to time constant of the impulse response of the filter. When we chop 
the filter impulse response at < igh value of T’, which requires a small value of M 
based on the relation M and 7” for a constant bit rate (M = T/T"), we conclude 
that more power is required to achieve the P, = 10-!* than chopping the impulse 
response at a lower value of T’. The reason for this is that the noise variance of 
the Fabry-Perot filter is inversely proportional to the value of M. As M decreases, 


more power penalty is paid to obtain the same bit error probability. 


When we compare the system employing a finite-time integrator and the sys- 


tem employing a Fabry-Perot filter, we conclude that the finite-time integrator has 
better performance with a high gain amplifier under low thermal noise values than 
the Fabry-Perot filter when we chop the impulse response of the filter at a high 
value of T’. For a larger value of M, we conclude that the Fabry-Perot filter has 
better performance with a high gain amplifier and low thermal noise values than 
the finite-time integrator. 

The system employing an optical preamplifier and a Fabry-Perot filter in a 
multi-channel environment suffers from intersymbol interference and adjacent chan- 


nel interference. When the channel spacing is large, the effect of ACI is small. 
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APPENDIX A 


Derivation of the Expected Value X; = E{X;} [33] 


Taking the expected of (4) we obtain 


x. - Eils(P} _ 1 e| 


: GP T°GP 


ies r,(r)dr 


t-T' 


(A.1) 


Substituting (1) into (Al) we have 
| 


= 1 
ie rasp | 
hf(G -1) 


= * le _B {itd} dry dt» y+ Seed (A.2) 


The laser phase noise 6(t) is characterized by a Wiener process [16] such that its 


T . 
as [VGP be) 4 ne(T) + jn,(r)| dr 


derivative is a zero mean white Gaussian process of PSD 278 where # is the laser 


linewidth. The variance of 0(t) is 27 ft and it can be shown that 
E {eil®(m)- ral = e774ln-7l (A.3) 


Substituting (A.3) into (A.2) yields 


i* (AT)? TGP ~~ 
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APPENDIX B 


Derivation of X; and G; 


X; is the squared envelope signal detected by photodetector with responsivity R 


and scaled by 1/T.RGP. It is defined as in (4) 
= 1 x : 2 
Xi = aaap f R|S,(t)(Pdt 
S(t) = Sa(t)=VEP bo f "a(t — r)dr 
h(t) = ce“, t>0 
S(t) = VGPbo(1-e*t), O<t<T 
7 1] T - 2 
X; = raaz | R|VGP boo (1 -e )| dt 


Xi = boo (1 = = (1 - eT) + = (1 - on) 


(B.1) 
(B.2) 
(B.3) 
(B.4) 
(B.5) 


(B.6) 


It is described in the first section, the noise n,(t) and n,(t) are independent over 


disjoint intervals and the random variables N., and N,, are independent identically 


distributed Gaussian random variables with zero mean. Their variances in (46) can 


be computed as follows, 


N-.(t) 


ae f i herintede 


A(t) = hoo(t)pri(t) 


1 O<t<T' 
pr(t) = 
Q = otherwise 
1 t 
Nt) = aes J beolt — 7)ne(r)pr(t — r)dr 


oF [ o halt=o nde 
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(B.7) 


(B.8) 


(B.9) 


1 nT’ 
N(2T’!) = —== h.o(nT’ — r)n.(1r)d B.10 
(nD) = ers [as bln! — P)nel rar (B.10) 
n(t) = = (B.11) 
of = of oa h2(nT’ — 1)dr (B.12) 
Ne Ne 9MGP Jin-yT! © 
using 7’ = nT’ — 7, 
N , , 
oh = oh, = stp [AL (r-dr") 
_ No [12 No [* 2 
= sap h Miner < sep ff Helen 
No 
on. =n, = IMGP© (B.13) 
where 
T! 
3 | h2,(r)dr (B.14) 


If we use the approximated RC filter of which impulse response h(t) given as 


in (B.3) 


oh, = OR, = sis (B.15) 
where 
es 
E = HE h?,(t)dt 
= [F (cers yat 
i] 
E = s(l-e-*) (B.16) 


We can approximate € as c/2 since the exponential term goes to zero due to the 
high value of 2cT”. 
(B.17) 
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Finally, the variable is found as 


: 2MGP 2 
Naphs(G — 1)e 
Dc ap''f 
ON. = ON, aMGP (B.18) 
We also know that this variance value can be represented by 
on, = oN, = 07 in pdf (54) (B.19) 
where N,, (the spontaneous emission factor) is unity and 
c= atl (B.20) 
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APPENDIX C 


Data signal in channel k: 


0 
b(t) = S$ by ce?’ pr(t — (T) (C.1) 
t=-L 
Data signal in channel 0: 
0 
bo(t) = 32 boupr(t —é7) (C.2) 
i=-Ly 


Bit in channel k in L" time interval (éT, (+ 1)T) 
bye € {0, elt} 
Bit in channel ¢ in i** time interval (7, (¢ + 1)T) 
bo; € {0,1} 
Detected bit in (0, T) 
boo € {0,1} 
depending on where bit zero or bit one is present. 


Desired channel = channel 0 


Adjacent channel = channel k, where 


k= ats --,-l, jee even integer 
2 2 
w, = Frequency spacing between channel k and channel 0 
QrkI 
Wr = T (C.3) 
¢, = Phase offset between channel k and channel 0 uniformly distributed over (0, 27) 
1 O<t<T 
pr(t) = 
0 otherwise 
(C.4) 
1 if<st<(@4+0T 
pr(t—tT’) = 
0 otherwise 
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The received signal at the output of the preamplifier and at the input of the Fabry- 


Perot filter (for channel ¢) 


M/ 
r(t) = VGP bo(t) + > VGP b;(t) (C.5) 
ku~M/2 


ky 
where P is the optical power and G is the optical amplifier gain. 


The output of the Fabry-Perot filter is 
ro(t)o ft h(t — r)r(r)dr (C.6) 
where A(t) is the equivalent lowpass impulse response of the Fabry-Perot filter of 


channel 0, 
h(t)=ce™* t>0 (€.7) 


= VGP [- h(t — r)bo(r)dr + VGP > fom h(t — 7)bg(7)dr 


ku —M/2 
k#0 


VGP too f h(t — r)pr(t)dr 


+VGP > bo ie h(t ~ r)pp(r — iT )dr 


t=-—Lo 


0 
+VGP x > bie a h(t — 1)" pr(r — €T)dr (C.8) 
ku —M/2 feb 
ko 


Since we are interested only in the detection interval 0 < t < T, we need to evaluate 


S(t) =ro(t), O<t<T 


M/2 
S(t) = VGP boo | h( h(t—r)dr + VOP bo. | h(t-—r)dr+VGP YS 


i=—Lo kn —M/2 


(i+1)T 


kgo 


¢41)T t 
b T)eI Vk" b h(t — r)e**" dr C.9 
(eB wef h(t —r)e | xo (t ) (C.9) 
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S(t) = Sp(t) + Sici(t) + Sacr(t) 


where 


Sa(t) 


Il 


VGP to [ "hi a)aF 


Stisr(t) = VGP > boa [on a(t — rr 
t=-Lo s 


M/2 -1 (¢+1)T . 
Saci(t) VGP > (| bre h(t — r)e?“*" dr 
t=-L 


fT 


kuM/2 
kyo 


t ; 
+t if h(t — rear] 


The evaluation of the integrals are as follows 


Salt) = VGP boo(1 — e-%), 0<t<T 
-1 
Sisi(t) = VGP 3S bo, (clttner _ gier) 
s=—Lo 
JGPe-t SY 
Saci(t) = VGP ce 3 iw 
kw—M/2 
k#0 


=-L 


Worst Cast Analysis 
bp, = b_, Ly = 00 


-1 
(| > bes ( elctiwe(C+1)T _ e(etsee) 
u 


0<t<T 


(C.10) 


(C.11) 


(C.12) 


(C.13) 


(C.14) 


(C.15) 


| + bio (elerivndt a ») 


S¥5(t)= V/GPbe-* O<t<T 
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(C.16) 


(C.17) 


bee = bio = b, L=o 


jwet 
S¥,,(t) = VGPb > = sa | OOSEST (C.18) 
ku —M/2 1+— 
c 
kyo 
with J = integer > 0 and wu, = 2xkI/T, 
x=R| ” 18;(t)/2at (C.19) 
ofS ; 


By using (C.14), (C.17), and (C.18), the evaluation of integrals gives the value of 
X as 


X = SREP [rreP (i {1-0-7 +5 spl — e*)] 


+82 [ (1 —e-*)] + bob. [a(1 26°? + e-*)| 


+ |r + 2 —(1 — e7°T)Re{b}(b_ — boo) 


/ 
> —ix))| (C20) 
ku —M/2 1+ (“) 


kyo 
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